Direct Detection of Kaluza-Klein Particles in Neutrino Telescopes 
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In theories with universal extra dimensions (UEDs), all standard model fields propagate in the 
bulk and the lightest state of the first Kaluza-Klein (KK) level can be made stable by imposing a 
Z-2 parity. We consider a framework where the lightest KK particle (LKP) is a neutral, extremely 
weakly interacting particle such as the first KK excitation of the graviton, while the next-to-lightest 
KK particle (NLKP) is the first KK mode of a charged right-handed lepton. In such a scenario, 
due to its very small couplings to the LKP, the NLKP is long-lived. We investigate the production 
of these particles from the interaction of high energy neutrinos with nucleons in the Earth, and 
determine the rate of NLKP events in neutrino telescopes. Using the Waxman-Bahcall limit for the 
neutrino flux, we find that the rate can be as large as a few hundreds events a year for realistic 
values of the NLKP mass. 

PACS numbers: 11.30.pb, 13.15+g, 12.60.jv, 95.30.Cq 



I. INTRODUCTION 

Although the standard model (SM) is a successful de- 
scription of the energy scales experimentally probed so 
far, it is expected that new physics will appear at the 
TeV scale. This is precisely the energy regime soon to be 
studied at the Large Hadron Collider (LHC). It is also the 
natural scale for the dynamical origin of electroweak sym- 
metry breaking, as well as for the solution of the hierar- 
chy problem. Typical solutions of these problems involve 
either symmetries (e.g. supersymmetry) , some dynam- 
ical mechanism of electroweak symmetry breaking (e.g. 
technicolor) , or a combination of symmetry and dynamics 
(e.g. little Higgs). In a somewhat different class are ex- 
tensions of the SM involving compact extra dimensions. 
In Large Extra Dimensions [l| only gravity propagates 
in the extra dimensional bulk, and the true fundamental 
scale of gravitation is 0(1) TeV. On the other hand, in 
theories with one curved extra dimension Q, gravity is 
weak at the TeV scale due to the warping produced by 
the bulk metric. 

Here we consider a more generic brand of extra di- 
mensional theories, universal extra dimensions (UEDs), 
where all fields propagate in the extra dimensional 
bulk 0. Its main motivation is phenomenological: if 
compact extra dimensions exist and all fields propagate 
in them, the inverse compactification radius could be 
just above the weak scale, setting the stage for a lot of 
new physics possibilities at the TeV scale. Furthermore, 
adding a mild assumption, the presence of a reflection 
symmetry leading to a ^-parity, UED theories are en- 
dowed with a candidate for dark matter: the lightest KK 
particle or LKP. 

Although at leading order the spectrum of each KK 
level is degenerate, it splits under radiative corrections, 
as well as when generic higher dimensional operators are 
taken into account In theories with one extra di- 



mension, if only the loop contributions coming from the 
physics below the cutoff are considered, one obtains the 
spectrum of the minimal UED standard model (mUED) 
of Ref. 4]. In this case, the LKP is most likely to be 
the first KK mode of the photon 7W. Other possibili- 
ties for the LKP include the first KK mode of the gravi- 
ton (?W and (in theories where neutrino masses are 
Dirac) the first KK excitation of the right-handed neu- 
trino Q. Other light particles include the KK ex- 
citation of a right-handed charged lepton, v- x \ and the 
charged Higgs KK mode Q. The splitting between the 
LKP and £^ is typically only a few GeV, depending on 
the choice of parameters [26|. 

The mUED spectrum is merely illustrative, and ul- 
traviolet physics contributions to boundary terms could 
significantly alter it, making for instance, v- 1 ' the NLKP, 
while either or remains the LKP. In such a 

scenario, the decay of the NLKP to the LKP would be 
highly suppressed, making the NLKP lifetime very large. 
We will consider this possibility in this paper. This sit- 
uation is analogous to what happens in some supersym- 
metric scenarios (e.g. gauge mediation) where the grav- 
itino is the lightest supersymmetric particle and a right 
handed charged slepton is the next to lightest one. The 
phenomenology associated to a long-lived £^ includes 
highly ionizing tracks at colliders. It also implies that 
PS 1 ^ can be produced by the interactions of high energy 
neutrinos with the earth and propagate through it until 
reaching a detector, in very close analogy to the case of 
NLSP sleptons studied in Refs. [I, E3, Cl E El • 

We will show that interactions of high energy neutri- 
nos (E„ > 10 5 GeV) with nucleons in the Earth will pro- 
duce pairs of NLKPs. The rate of events will allow the 
discovery of the latter in km 3 neutrino telescopes. This 
analysis follows the same steps as for NLSP detection 
[i| [l2j]- The crucial observation is the same as for the 
NLSP, the small NLKP production cross section is com- 
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pensated by its large range. The NLKP loses much less 
energy while traveling through the Earth when compared 
to SM leptons. This allows the detection of NLKPs that 
are produced far away from the detector. 

As the NLKPs are produced in pairs, the main back- 
ground consists of di-muon events. We will show that 
there are at least two ways to separate these from the 
signal. For lower mass NLKPs, the measured energy 
spectrum will have a bump in the region from 10 3 to 
10 4 GeV due to the fact that the energy loss in the de- 
tector will resemble the one from lower energy muons. In 
addition, for both low or high mass NLKPs, the separa- 
tion between the pair that crosses the detector will be 
larger than the one for the di-muon pair, and will allow 
to distinguish the signal from the background. 

This paper is organized as follows: we first determine 
the NLKP production cross section; in Section ITU we de- 
scribe the NLKP energy loss while traveling through the 
Earth; the analysis of the signal and comparison with the 
background are discussed in Section HVl and the conclu- 
sions follow in the last section. 



II. NLKP PRODUCTION 

In this section, we compute the production cross sec- 
tion for the NLKP pair. Due to the presence of the im- 
parity, all KK modes produced will eventually cascade 
down to a NLKP. Since KK modes are produced in pairs, 
KK production initiated by vN scattering will result in a 
pair of NLKPs. The dominant process for v — TV-initiated 
KK production involves the t-channel production of a left 
handed lepton KK mode L\ (with generation index i) 
and a quark KK mode (Q^ ) via a gauge boson KK mode 
(W {1) ). This process is analogous to the charged current 
(CC) in the SM. We also include the subdominant pro- 
cess which is analogous to the neutral current process in 
the SM. This involves the exchange of a neutral gauge 
boson KK mode These processes are shown in 

Figure [T] 

The neutrino, which is always left-handed, can inter- 
act with a left handed down- type quark (a) or with a 
right-handed up- type antiquark (b). This results in the 
partonic cross sections: 

dffW W (^-<(D -^q(d) 

dt ~ sm 4 6 w s 2 (i 2 -M^ (1) ) 2 ( ' ' 

efo-O) 87ra 2 

—r- = — 3 x ( 2 - 2 ) 

dt sin 9w 

m l(i) m Q(D +u 2 -u (m 2 (1) + m 2 QW ) 

a 2 (t 2 -M^ (1) ) 2 

(2.3) 

where s, t and u are the usual Mandelstam variables, 
and m L (i), mqw and M w w are the Lf \ and 




(c) Q U (d) Q 

FIG. 1: Feynman diagrams for KK mode production in uN 
collisions. Charged current (charged gauge boson KK mode) 
interactions: (a) Left-left interaction requiring the insertion of 
the gauge KK mode mass in the t-channel line, (b) Left-right 
interaction. Neutral current: (c), (d). There are analogous 
diagrams for anti-neutrinos as well as for strange and charm 
initial quarks. 

masses, repectively. The subdominant neutral gauge bo- 
son KK mode (Z^ 1 ') exchange is shown in Figure [1] (c)- 

(d). Each of these processes will produce a and a 
Q^ 1 ) and both of these particles will promptly produce a 
decay chain ending with a 

Bounds from direct searches from the Tevatron, as 
well as from electroweak precision constrains Q , result in 
R- 1 > 300 GeV for 5D, while for 6D is R- 1 > 500 GeV. 
We will assume three illustrative values for the NLKP 
mass: 300, 600 and 900 GeV. Finally, we need to specify 
the cutoff of the theory. Using naive dimensional analy- 
sis, we find for the 5D case, that AR ~ 8ir/g 2 , where g 
is the strong coupling constant. We then take AR ~ 20 
as a sharp cutoff, i.e. we include no contributions from 
energies above 20/ R. In order to evaluate the uncer- 
tainty introduced by this procedure in the cross section 
calculation, we scanned values of AR up to 30, with no 
significant effects in the results. 

The NLKP production cross section is shown in Fig- 
ure 2 as a function of the neutrino energy. For com- 
parison, the SM charged current (top gray curve) and 
the di-muon (solid red curve) background cross sections 
are also shown. As expected, the NLKP production cross 
sections (ctnlkp) are significantly lower than the SM one. 
However, depending on the neutrino energy and the 
mass, it can be larger than the di-muon background. In 
the next section we will show that the the fact that ctnlkp 
is rather suppressed (as compared to the SM one) will be 
compensated by the sizable NLKP range resulting from 
the combination of its long lifetime and small energy loss. 

It is also interesting to compare the ctnlkp to the NLSP 
production as obtained in Ref. The NLKP produc- 
tion is significantly larger than the one for NLSPs, which 
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FIG. 2: vN cross sections vs. the energy of the incident neu- 
trino. The violet dashed, blue circled and black crossed lines 
correspond respectively to 300, 600 and 900 GeV NLKPs. 
The top gray curve corresponds to the SM charged current 
interactions and the full red one to the di-muon background. 



translates into a larger number of events at the detector, 
as we will see below. 



In order to determine the NLKP energy losses, we fol- 
low closely the calculations done for NLSP propagation 
in Refs. [l2j and [l7|. Radiation losses dominate above 
a propagating energy of 1 TeV. Among them, pair pro- 
duction and bremsstrahlung are less important for the 
NLKP when compared to photo-nuclear interactions, as 
can be seen in Figure 3. Even so, the energy loss due 
to photo-nuclear interactions is suppressed by the NLKP 
mass. As mentioned in Ref. @ and shown explicitely in 
Refs. [12] and [17], the important energy region for this 
process is the one at low photon virtuality Q 2 . The rea- 
son is that the structure function involved in the process 
is determined by a cross section which is dominated by 
physics at low Q 2 ~ 1 GeV 2 . However, due to the large 
NLKP mass, the minimum value for the photon virtuality 
will be larger, therefore avoiding the effects of resonances 
and other nonperturbative processes which occur at lower 
Q 2 . This is in contrast to the case of the r lepton, where 
the resonant region still dominates and results in a much 
larger photo-nuclear energy loss. 

Figure 3 shows the radiation loss term of eq. (|3.4|) ver- 
sus neutrino energy for muons, taus and the 300 GeV 
NLKP. As expected, the photo-nuclear process domi- 
nates the NLKP radiation loss. However, it is still quite 
suppressed due to the NLKP heavy mass and the total 
energy loss is still considerably below the one for SM lep- 
tons. Energy suppression will be enhanced for heavier 
NLKP mass. 

We then conclude that the NLKP energy loss is quite 
suppressed in comparison with SM leptons. As we will 
see below, this means that its range through the Earth 
is much larger, allowing for the detection of NLKPs that 
have been produced hundreds or even thousands of kilo- 
meters from the detector. 



III. NLKP ENERGY LOSS 

After production, the NLKPs lose energy due to ion- 
ization and radiation processes. The average energy loss 
is given by [T3 |: 

— = a(E) + b(E) E , (3.4) 

ax 

where a(E) represents the ionization losses, and b(E) the 
contributions from different radiation processes. The lat- 
ter includes bremsstrahlung, pair production and photo- 
nuclear interactions. There is also energy loss due to 
weak interactions, but this will only be of importance 
at very high energy [LSI ], and we will neglect it for the 
remaining of this work. 

At the high energies where NLKPs can be produced, 
radiation losses dominate over ionization. Among radia- 
tion processes, both pair production and bremsstrahlung 
become less important for heavy particles. Although 
photonuclear processes dominate tau lepton propagation 
losses [lfl [HI , a mass suppression will occur for leptons 
of much heavier masses I12L Il7ll . 



IV. NLKP SIGNALS AND RATE IN NEUTRINO 
TELESCOPES 

A. Neutrino Flux 

The NLKP event rate in neutrino telescopes depends 
on the incoming neutrino flux. This is largely determined 
by the high energy cosmic ray spectrum [19j . There are 
other potentially relevant sources of the neutrino flux, 
such as atmospheric charm production (20| . For the pur- 
pose of this work we will neglect these other contribu- 
tions, only considering the flux of cosmic neutrinos, for 
which we use two alternative estimates: the work of Wax- 
man and Bahcall (WB) [2l| and the one of Manhcim, 
Proterhoe and Rachen (MPR) [22j • The integrated num- 
ber of events resulting from the MPR limit is consider- 
ably larger than the WB. We find our NLKP rates assum- 
ing each of these limits as our incoming neutrino flux. All 
plots are produced assuming the WB limit as our neu- 
trino flux. 

Waxman and Bahcall fix the cosmic ray spectrum to 
a power law curve with spectral index —2. The neutrino 
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FIG. 3: Radiation energy loss b(E) parameter due to pair 
production, bremsstrahlung and photonuclear processes for 
muon, tau and a 300 GeV NLKP. The plot labeled "Tot" 
is the sum of all contributions. The muon, tau and NLKP 
curves are as labeled in the first plot. Heavier NLKPs will 
have lower b(E) parameters. 



FIG. 4: Upper bound on differential neutrino flux as calcu- 
lated by WB [U (shaded area) and MPR ^22] (red lines). The 
WB limit ranges from the limit with no cosmological evolution 
(upper edge) and with cosmological evolution (lower edge). 
The MPR limit is shown for optically thin sources (red dot- 
ted line) as used in this paper and for optically thick sources 
(red dot-dashed line). 



upper limit is deduced assuming that each nucleon will 
interact with photons and produce a pion. The charged 
pions will then decay producing neutrinos. Their argu- 
ment requires that the sources are optically thin, which 
means that most of the protons escape the source with- 
out interacting. As a result, the neutrino upper limit is 
given by 



d<j) v 
~dE 



WB 



(1 -4) x 10 



-s 



-GeV cm 2 s 1 sr- 



1 



where the range in the coefficient depends on the cosmo- 
logical evolution of the sources. The evolution accounts 
for the source activity and redshift energy loss due to the 
cosmological expansion. We take the upper end as the 
neutrino flux incoming through the Earth. 

On the other hand, instead of taking a fixed power law 
behaviour for all cosmic ray spectrum, Manheim, Proter- 
hoe and Rachen determine the spectrum at each energy 
directly from data. Here we consider the limit MPR ob- 
tain assuming optically thin sources, although they also 
determine a limit for optically thick sources (See com- 
ments about optically thick sources in [Hj]). Figure 4 
shows both WB and MPR limits for the muon plus anti- 
muon neutrino flux. 

As seen in Section [Til the NLKP production is inde- 



pendent of the initial neutrino flavor. For this reason we 
consider both electron and muon neutrinos, and neutrino 
mixing does not affect our results. 



B. NLKP Signals 

We now have all the ingredients to determine the 
NLKP rate at neutrino telescopes. In order to under- 
stand the signal in detail, we performed a Monte Carlo 
simulation generating approximately 30,000 events for 
each NLKP mass (300, 600 and 900 GeV). 

Once the incoming neutrino flux is determined, an in- 
teraction point is randomly chosen based on the NLKP 
production probability. This results from a convolution 
of the neutrino survival probability with the probability 
of interacting and producing a NLKP. The neutrino sur- 
vival probability P$ is given by exp(J ndl), where n is the 
Earth number density and I is the distance the neutrino 
travels. We use the Earth density profile as described in 

[nm. 

The primary particles {L^ and Q^ 1 ' produced in the 
neutrino interaction) angular distribution at the CM is 
randomly determined based on the differential produc- 
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FIG. 5: NLKP pair event energy distribution per km 2 , per 
year, at the detector. Plus violet line corresponds to 300 
GeV; blue circled line to 600 GeV and crossed black line to 
900 GeV NLKP. For reference the neutrino flux at earth (full 
black line); and the /x (dotted green line) and the di-muon 
(squared red line) flux through the detector are also shown. 
In all cases we make use of the WB limit for the neutrino flux. 



tion cross section. The center of mass (CM) angular dis- 
tribution of the two NLKPs produced is assumed to be 
the same as the one between the two primary particles. 
This is a good approximation [i~2| for events with energy 
well above the production threshold where most of the 
event rate comes from. The events close to the produc- 
tion threshold have a broader angular distribution. These 
would enhance the separation differences between signal 
and background and therefore make our results conser- 
vative. The CM angular distribution is then boosted to 
the laboratory frame. 

Once the NLKPs are produced their propagation 
through the Earth is simulated. Their energy loss - which 
is mass dependent - is taken into account. The NLKP 
energy distribution as a function of neutrino energy is 
shown in Figured] As can be seen, the 300 and 600 GeV 
NLKP event rate is much larger than the muon's for en- 
ergies above NLKP production threshold. The 900 GeV 
NLKP rate will be comparable with that for muons, but 
still larger than the di-muon background rate. 

Although these are rather large rates, they do not di- 
rectly translate into observed NLKPs due to the fact that 
NLKPs are hard to identify. Neutrino telescopes measure 
their energy in two ways [TH, (25[: low energy events (be- 



TABLE I: Number of events per km 2 per year for different 
NLKP masses and neutrino fluxes at the Earth. The NLKP 
masses are 300, 600 and 900 GeV. The number of NLKP 
events are given for energies above threshold for production 
of a and a Q' 1 ' while the muon rate for energies above 
1000 GeV. The column corresponds to the di-muon 

background. No cuts were applied at this stage. 









(300) 


r(i)T (1) 
(600) 


(900) 


WB 
MPR 


552 
39654 


30 
1914 


489 
1476 


21 

47 


3 
5 



low ~ 100 GeV) have their energy reconstructed from 
the track length, whereas for the more energetic ones 
the energy is reconstructed from the amount of Cerenkov 
light deposited in the photomultiplier tubes. Taking the 
Cerenkov radiation as proportional to the amount of de- 
posited energy in the detector is a good approximation 
for SM leptons. But the NLKPs lose a lot less energy 
than SM leptons. Thus, if a NLKP track is assumed to 
be a SM lepton such as a muon, it will be assigned a 
much lower energy as such. For this reason and in order 
to compare event rates, the muon rate must be integrated 
from energies lower than the KK production threshold. 

Table [T] shows the event rate per year per km 2 both for 
the WB flux, as well as for the MPR optically thin flux. 
The numbers are clearly encouraging for km 3 neutrino 
telescopes. Two features will be important to distiguish 
the signal from the background : the separation between 
the pair of particles inside the detector; and - for lower 
mass NLKPs - a bump in the energy spectrum will ap- 
pear. These features will be discussed at the end of this 
section. 



1. Di-muon Background 

Due to their large boost most NLKP pairs go through 
the detector in two well separated and approximately 
parallel tracks. Events well separated are produced far 
from the detector and as the production angle between 
them is small the tracks will be almost parallel. There- 
fore, single muons can be eliminated by a two track 
requirement. The main remaining background are di- 
muons. These are produced from charmed hadrons from 
the following chain : 

vN — > fi~ H c — > /j," /i + H x v , 

where H c is a charm hadron produced from a muon neu- 
trino CC interaction and H x can be either a strange or 
non-strange hadron. 

The cross section for charm production from a neu- 
trino interaction was calculated in Ref. as well as 
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the di-muon energy loss, propagation and separation at 
the detector. In what follows we reproduce these results, 
and compare with the NLKP signal. 



2. NLKPs Separation 

The separation between the NLKPs will be given by 
the distance traveled times the angle (9) between the pair 
in the laboratory frame. As the boost from CM to lab 
is large, 9 is very small. However, this is compensated 
by the production being far away from the detector. The 
production point being typically a few 1000 km from the 
detector and 9 ~ 10~ 4 — the separation between the 
two NLKPs will be a few tens to a few hundred meters. 

On the other hand, di-muon events have to be pro- 
duced close to the detector, otherwise they lose all their 
energy before arriving at it. For this reason their sepa- 
ration is typically smaller than the one for most of the 
signal events. 

The separation distribution for each NLKP mass at 
the detector is shown in Figure [6] The simulated detec- 
tor is placed at the same depth as the IceCube telescope 
[25l | . We also show the di-muon background separation 
for comparison. While the dimuon separation is at most 
~ 100 m, the pair of NLKP can be more than 100 me- 
ters apart. For instance, for a 300 GeV NLKP, 52% of 
the events are more than 50 m apart and 28% are more 
than 100 m apart. The di-muon background has only 8% 
with more than 50 m and 1.3% with more than 100 m 
separation. The 600 and 900 GeV NLKPs have both 
around 60% of events with more than 50 m separation 
and around 42% with more than 100 m separation. 

In order to estimate the statistical significance of the 
separation cut, we determine the Sj\f~B ratio, where S 
and B are respectively the number of signal and back- 
ground events. We find that for the 300 GeV NLKP, a 
requirement that the pair of NLKPs are at least 10 me- 
ters apart will yield a significance of 85, ie, 436 of the 489 
NLKPs will be more than 10 meters apart, while only 25 
di-muons will have more than 10 meters separation. For 
the 600 GeV NLKP, a requirement of 86 m separation 
will allow a 5a significance in one year, with 9 signal 
events and 3 di-muons. For the 900 GeV, the separation 
is harder, a 2a significance can be achieved in a year with 
a separation cut of 150 m, while a 5a significance needs 
5 years to be achieved. 



3. The NLKP Bump 

Another feature of the NLKP signal comes from the 
fact that these particles lose less energy than a SM lep- 
ton. This implies that NLKPs will have their energy 
reconstructed as if they where lower energy leptons. Fig- 
ure [7] shows both NLKP and di-muon simulated energy 
distribution as they arrive at the detector. Although 
the NLKPs are more energetic than the di-muons, the 
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FIG. 6: Track separation distribution between NLKP pair 
(top) and for the di-muon background (bottom). 



energy deposited in the PMTs will resemble lower en- 
ergy muons and therefore they have to be compared with 
them. However this will generate a sizeable excess in the 
reconstructed energy spectrum, at least if the number of 
NLKP events is large enough. In order to understand 
how this feature will change the reconstructed energy 
spectrum, we simulate the reconstructed energy by tak- 
ing all NLKPs as muons. This was done by determining 
the deposited energy in the detector and reconstructing 
this energy as if deposited by a muon. These events were 
then added to the SM muon energy spectrum. The con- 
sequence is that the high energy NLKPs will be recon- 
structed as lower energy events that will end up as a 
bump around energies of TeVs. 

Figure[8]shows the energy distribution of the muon flux 
through the detector (top plot, blue circles) and the same 
distribution with the addition of 300 GeV NLKPs recon- 
structed as muons. A visible "crown" with few events in 
each energy bin in the 1 to 100 TeV region will clearly 
indicate the presence of KK particles. This feature will 
be enhanced when the NLKPs are included in the di- 
muon energy spectrum (bottom plot). When the signal 
is reconstructed as di-muons, a pronounced crown shows 
up in the reconstructed energy spectrum. This feature is 
observable for NLKPs in the lower mass range, since the 
rate of higher mass NLKPs would not be large enough 
to observably enhance the spectrum in the lower energy 
region. We expect this feature to be observable up to 
NLKP masses of about ~ 600 GeV. Thus, for these lower 
mass NLKPs there will be two distinct ways to separate 
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FIG. 7: Arrival energy distribution of the NLKP at the de- 
tector for m (i) = 300, 600 and 900 GeV. Also shown is the 

R 

arrival distribution for the di-muon background. 



the signal from the main background. 



We have made a detailed study of the background and 
the signal, and shown that the track separation of NLKPs 
is a good discriminant with respect to the di-muon back- 
ground. In addition, for the case of smaller NLKPs 
masses, we have shown that the NLKP signal results in 
a bump in the detected di-muon spectrum, since NLKPs 
lose energy similarly to lower energy muons. Combining 
this feature with the characteristic track separation of the 
signal tracks should enhance the statistical significance of 
a potential signal. 
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V. CONCLUSIONS 

We have shown that in a UED scenario where the 
NLKP is a the first KK mode of a right handed charged 
lepton, neutrino telescopes such as IceCube will be able 
to directly observe these i^'s up to masses of several 
hundred GeV, perhaps even 1 TeV. This complements 
hadron collider searches, where signals for this UED sce- 
nario would consist of large missing energy, and per- 
haps one or two highly-ionizing tracks. The similarity 
of the UED signals with the analogous supersymmetric 
scenario, for instance with gravitino dark matter and a 
slepton NLSP, can make the identification of the under- 
lying theory difficult. On the other hand, the event rate 
at neutrino telescopes coming from this UED scenario 
is considerably higher than the one resulting from the 
supersymmetric case and studied in Refs. 0, [l2j ■ 



FIG. 8: Top: Energy distribution of the muon flux through 
the detector (blue circles) and the same flux plus the 300 GeV 
NLKP spectrum reconstructed as muons (violet plus line). 
Bottom: Same as above but now using the di-muon flux 
through the detector. 
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